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During the first 500 million years of cosmic history, the first 
stars and galaxies formed and seeded the cosmos with heavy ele- 
ments [e.g., 1-4]. These early galaxies illuminated the transition 
from the cosmic “dark ages" to the reionization of the intergalactic 
medium. This transitional period has been largely inaccessible to 
direct observation until the recent commissioning of JWST, which 
has extended our observational reach into that epoch. Excitingly, 
the first JWST science observations uncovered a surprisingly high 
abundance of early star-forming galaxies [e.g., 5-12]. However, the 
distances (redshifts) of these galaxies were, by necessity, estimated 
from multi-band photometry. Photometric redshifts, while generally 
robust, can suffer from uncertainties and/or degeneracies. Spectro- 
scopic measurements of the precise redshifts are required to validate 
these sources and to reliably quantify their space densities, stellar 
masses, and star formation rates, which provide powerful constraints 
on galaxy formation models and cosmology [13]. Here we present the 
results of JWST follow-up spectroscopy of a small sample of galaxies 
suspected to be amongst the most distant yet observed. We confirm 
redshifts z > 10 for two galaxies, including one of the first bright 
JWST-discovered candidates with z — 11.4, and show that another 
galaxy with suggested z = 16 instead has z = 4.9, with strong emis- 
sion lines that mimic the expected colors of more distant objects. 
These results reinforce the evidence for the rapid production of lumi- 
nous galaxies in the very young Universe, while also highlighting the 
necessity of spectroscopic verification for remarkable candidates. 

The spectroscopic targets were selected from JWST/NIRCam [14-16] 
imaging obtained through the Cosmic Evolution Early Release Science 
(CEERS) Survey [17] because their colors at wavelengths of 1 — 5 microns 
matched predictions for galaxies with redshifts z > 9 [12], including the key 
signature of a “break” in flux density at wavelengths shorter than that of red- 
shifted hydrogen Lyman a (Lya; rest-frame wavelength Arest = 1216À) due to 
the largely-neutral intergalactic medium (IGM) expected at early times. Here 
we focus attention on two bright extreme-redshift candidates, one with a pho- 
tometric redshift z ~ 11.1 (J141946.36+525632.8, aka Maisie’s Galaxy [7], here 
MSA ID 1; see Methods for updated photometric redshift measurement) and 
another with a photometric redshift z ~ 16.4 (CEERS-93316 [8]; here MSA 
ID 0). At their estimated redshifts, these two galaxies would have UV lumi- 
nosities ~ 4 and ~ 15 times brighter than the most distant spectroscopically 
confirmed galaxy to date (z = 13.2 [18]). We show these photometric redshift 
estimates in Figure 1. 

Although independent analyses by different teams [7, 8, 12, 19] identified 
these two candidates using 7-band JWST/NIRCam photometry, it is possi- 
ble for foreground objects to infiltrate high-redshift samples due to unusual 
colors resulting from older stellar populations, strong nebular emission lines, 
and/or dust reddening [8, 20, 21]. A tentative detection at 850 um of dust 
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Fig. 1 The likelihood (bottom) and x? goodness-of-fit (top) distribution functions for the 
photometric redshift, derived from the JWST CEERS photometry. Results for Maisie’s 
Galaxy [7, 12] are in blue, and the z ~ 16.5 candidate galaxy discovered by [8] are in red. The 
line shading denotes the results using the photometric measurements from each team. The 
results between the two methods are broadly similar, though our use of HST photometry 
(see Methods) results in a redshift closer to the spectroscopic value due to a weak detection 
in the HST/WFC3 1 jum band. The spectroscopic measurements of z = 11.44 and z = 4.912 
from this paper are shown by the thick vertical lines. 


emission in the vicinity of CEERS-93316 [22] could support the dust reddening 
hypothesis, but we have investigated this using new observations at 1.1 mm 
with the NOEMA interferometer and localized the dust emission to another 
galaxy approximately 1.5 arcsec away (see Methods for details). Nonetheless, 
reddening could still affect the observed colors of CEERS-93316. 

We observed both galaxies with JWST/NIRSpec [23] using its multi- 
object spectroscopy mode and prism disperser, whose low spectral resolution 
at wavelengths shorter than 2.5 jum affords high sensitivity for detecting faint 
continuum emission and spectral breaks. Details of the observations and data 
analysis, as well as results from simultaneous observations of several other 
candidate z > 9 galaxies, are summarized in the Methods section. 

The spectrum of Maisie’s Galaxy (Figure 2) shows well-detected contin- 
uum emission with a sharp break, below which no significant flux is detected. 
The most likely interpretations of this feature are the Lya IGM break (Arest 
= 1216A) or the Balmer break (Arest = 3646À); the latter is prominent in 
galaxies where star formation ceased more than 100 Myr before the time of 
observation. Integrating the spectra in wavelength intervals redder and bluer 
than the observed breaks, we measure a 3c limit on the flux density ratio 
> 2.6, inconsistent with values expected for Balmer breaks [18] (see Methods). 
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Fig. 2 a) A color image of Maisie’s Galaxy, showing approximately the true rest-ultraviolet 
color (made from the F277W, F356W, and F444W images). The blue color is indicative 
of low dust attenuation and a fairly young stellar population. b) The region of Maisie’s 
Galaxy’s spectrum around the detected Lya break. The top panel shows the 2D spectrum as 
the bright trace, with the negative “shadows” representing the subtractions of the individual 
nods. The red horizontal lines show the spatial extent over which the 1D spectrum, shown 
in the bottom panel, was extracted. A clear spectral break is seen, below which there is no 
significant flux. We interpret this as the Lya break at a redshift of 11.44, which represents 
a time 390 Myr after the Big Bang. 


1 kpc 


We used several methods to measure the wavelength of the spectral break 
(see Methods and Table 2). Interpreting this feature as the Lya break (Fig. 2b) 
yields a redshift z = 11.44+9:03, consistent with photometric redshift esti- 
mates. Performing stellar population modeling (see Methods for details) on the 
observed photometry of this galaxy, fixing the redshift to this spectroscopic 
value, we confirm earlier results [7, 12] that Maisie's Galaxy is a modestly mas- 
sive (log;; M./Mg = 8.4*01), dust-poor (Ay = 0.07430? mag) galaxy, with a 
high star-formation rate for its stellar mass (log )(SFR/M.) = —8.0*03 yr- 1). 

NIRSpec observations of star-forming galaxies at z < 9.6 frequently detect 
strong [O 111] and HG emission lines [24-27] that shift beyond the instrument’s 
spectral coverage for redshifts z > 9.6. Emission lines at blue optical and 
ultraviolet rest-frame wavelengths are typically weaker and harder to detect 
in low signal-to-noise ratio (S/N) spectra. The observed spectrum of Maisie’s 
Galaxy shows no strong emission lines (a strong feature potentially consis- 
tent with emission from singly-ionized oxygen, [O11] rest-frame wavelength 
Arest = 3727A, observed at 4.6m, was discarded as a likely instrumental 
artifact — see the Methods section for further explanation). Exploring the 
rest-frame ultraviolet portion of the spectrum shows several weakly positive 
features. Specifically, there are positive emission features from doubly-ionized 
oxygen ([O rii], rest-frame wavelength Arest = 1666À) and doubly-ionized car- 
bon (Cin, Arest = 1909A) at a redshift of z = 11.383. However, we estimate 
that these emission lines are significant at only the 2-30 level. Thus we do not 
consider them robust detections, though they would corroborate the Lya-break 
derived redshift. 

The spectrum of CEERS-93316 (Figure 3) shows emission lines from the 
hydrogen Balmer series, [O i1], [O 11], and [S11] (Arest = 6718,6733A) at a red- 
shift z = 4.912+0.001. This is the same redshift as the galaxy CEERS-DSFG-1 
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Fig. 3 The NIRSpec spectrum of CEERS-93316. Similar to Figure 2, Panel a shows the 2D 
spectrum with the red lines denoting the spatial extraction limits, and Panel c shows the 1D 
extraction. Clear emission lines can be seen in both panels, which we label in the bottom 
panel. The bottom panel also illustrates the wavelengths covered by the JWST/NIRCam 
filter set used to select this galaxy, highlighting how [O 11] contamination in F277W (green) 
and Ho contamination in all of F356W (yellow), F410M (orange) and F444W (red) conspired 
to mimic the signature of both a very high-redshift Lya break and a blue continuum at 
longer wavelengths. This unusual (but educational) situation is detailed further in Figure 4. 
Panel b shows a F150W+F200W+F277W color image of this galaxy. 


(MSA ID 2), which lies 26 arcsec away and was previously discussed by [22] as 
a moderate-redshift galaxy with colors that could mimic a high-redshift Lya 
break. This redshift is also consistent with the alternative lower-redshift solu- 
tion proposed for CEERS-93316 [8, 20-22], and the same redshift as (at least) 
two other neighboring galaxies (see Methods). The continuum emission of 
CEERS-93316 fades at wavelengths shorter than 2.5 um (rest frame S 0.4 um), 
although flux is significantly detected in the spectrum down to A ~ 1.2 um. 
The strong Ha (Arest = 6564A) and [O 111] lines have large rest-frame equiva- 
lent widths (790 + 140A and 1020 + 180A respectively), and hence contribute 
significantly to the flux measured by NIRCam through the F277W, F356W, 
F410M, and F444W filters (specifically, at this precise redshift, Ha emission 
actually contributes to the flux in the last three of those filters; Figure 3). The 
combination of the strong lines and the red continuum leads to the red col- 
ors measured between filters F150W, F200W, and F277W and the blue colors 
from F277W to F444W (Figure 1b). 

The significant reddening implied by the relative lack of observed UV con- 
tinuum flux is corroborated by the large Balmer decrement measured from the 
spectrum (Ha/H8 = 5.6 + 0.6), which implies - despite the large equivalent 
widths of the lines - a substantial nebular color excess of E(B — V)nebular = 
0.71 + 0.12 assuming a Cardelli extinction law (see Methods). Indeed such 
scenarios have been considered as alternative explanations for the colors of 
CEERS-93316 [8, 20, 22]. An active galactic nucleus could also be responsible 
for strong line emission in an otherwise older galaxy [e.g., 28, 29]. 
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Fig. 4 Spectral energy distributions for Maisie's Galaxy (a), and for CEERS-93316 (b). 
Both panels show the observed photometry compared to the best-fitting stellar population 
model (black) and the observed spectra. In the bottom panel, the model is deliberately 
shown at z — 16.2 (consistent with the original photometric redshift) to demonstrate the 
way in which strong emission lines impact the observed photometry, conspiring in this spe- 
cific (and somewhat unusual) case to replicate the presence of the Lya break at z ~ 16 and 
an apparently blue continuum at longer wavelengths. The squares show the observed pho- 
tometry with the spectroscopic emission line contributions (shaded brown in the spectrum) 


removed. 


'The example of CEERS-93316 is a reminder that the prevalence of very 
strong nebular line emission in early galaxies can influence broad-band pho- 
tometry with instruments such as JWST/NIRCam and mimic the colors 
predicted for much more distant objects [22]. Exceptional candidates for galax- 
ies at extreme redshifts require spectroscopic confirmation. However, we note 
that this is a fairly extreme case, as it is only at very specific redshifts that 
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Table 1: Summary of Galaxy Properties 


Property Maisie's Galaxy CEERS-93316 
RA (J2000) 214.943152 214.914550 
Dec (J2000) 52.942442 52.943023 
Photometric Redshift ^ 11.08*035 16.451013 
Spectroscopic Redshift 11.44* 009 4.912 + 0.001 
TBigBang (Myr) 391 1195 

log( M. / Mc) gus 8.5401 

Ay (mag) DT E 1460729 
log(sSFR/yr-!) -8.0*02 —8.0*01 


Table 1 All uncertainties listed are 68% confidence intervals. The photometric redshift 
calculation is described in the Methods and has been updated for this work. We list here 
the spectroscopic redshift derived from the method with the smallest uncertainty for a 
given source; see Methods for a full discussion. 7pigpang denotes the age of the universe at 
the spectroscopic redshift. The stellar population properties are derived from Bagpipes. 


strong emission lines can fully mimic a Lya break with multiple broad-band 
filter detections. 

Fortunately JWST/NIRSpec is a powerful tool for spectroscopic analysis 
and redshift measurement, as demonstrated by our observations of Maisie’s 
Galaxy and three other galaxies at 7.9 < z < 10.1 (discussed in Meth- 
ods). NIRSpec observations have now confirmed redshifts for 9 galaxies at 
z 2 10 [18, 30, Arrabal Haro in preparation], including two with Muy < —20, 
unexpectedly bright compared to most pre-JWST predictions. 

These confirmations are crucial to improving our physical understanding 
of the first phases of galaxy formation. Early analyses of candidate galaxies at 
z = 10-12 have shown that their observed abundance is significantly higher 
than predicted by nearly all theoretical models [e.g., 7, 10, 12]. This could 
indicate that specific assumptions in these models need to be modified at very 
high redshifts, including an increased efficiency of converting gas into stars 
at the earliest times or even a higher characteristic star mass in these early 
epochs. Such a “top-heavy” initial mass function has indeed been predicted to 
be present in early galaxies, driven by a combination of their chemically simple 
makeup [e.g., 31, 32] as well as a rising cosmic microwave background tem- 
perature floor [e.g., 33]. However, before any such unusual physical properties 
can be robustly inferred for early galaxies, their photometric redshift estimates 
must be verified. The NIRSpec confirmation of Maisie's Galaxy at z — 11.4 
reported here, combined with growing numbers of z = 10 spectroscopic con- 
firmations in the literature [e.g., 18], show that well-established methods of 
photometrically selecting galaxies work robustly at z ~ 10-12. 


Methods 


Near-Infrared Imaging data and Source Selection 


We select high redshift galaxy candidates from the first epoch of imaging 
from the Cosmic Evolution Early Release Science (CEERS; PI: S. Finkelstein 
[7]) survey. These data consist of four roughly contiguous pointings with the 
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NIRCam instrument on JWST, covering an area of ~35 arcmin?. The full area 
is covered by six broad-band filters (F115W, F150W, F200W, F277W, F356W, 
F444W) and one medium-band filter (F410M). We use the v0.5 data release 
from the CEERS team! [17]. These reductions use the STScI JWST pipeline 
with several custom modifications, correcting for various issues, including 1/ f 
noise, wisps, and snowballs, with updated photometric calibration from late 
2022. The 50 photometric depths of these data for compact source detection 
are ~29-29.2 AB magnitudes in F115W, F150W, F200W, F277W and F356W, 
and ~0.7 and ~0.5 mag shallower in F410M and F444W, respectively. 

Spectroscopic observations targeted three primary galaxy candidates: 
CEERS J141946.36+525632.8, aka “Maisie’s Galaxy” [7], with photometric 
redshift z ~ 11.1, “CEERS-93316” with photometric redshift z ~ 16.4 [8], and 
CEERS-DSFG-1, a JWST-detected galaxy with millimeter-wavelength emis- 
sion that is likely to arise from dust-obscured star formation at z ~ 5, but 
whose near-infrared colors could be interpreted as suggesting z > 10 [22]. 
'These object names come from the original papers reporting these objects 
rather than from a homogeneous source catalog. Essential parameters of the 
first two galaxies, including results from the present analysis, are summarized 
in Table 1. Detailed analysis of CEERS-DSFG-1 will be presented elsewhere. 

Photometry was optimized for accurate measures of colors, generally fol- 
lowing [12], with a few key updates to improve measured colors, particularly 
between HST and JWST. Specifically, all bands with PSFs smaller than that of 
F277W (i.e., ACS F606W, F814W, and NIRCam F115W, F150W, and F200W) 
were convolved to match the F277W PSF. For the remaining images, fluxes 
were measured in the native images, but a correction was applied on a per- 
source basis as the ratio of the flux in the native F277W image to that of the 
F277W image convolved to match the PSF in a given larger-PSF image. We 
then used simulation-based corrections to derive accurate total fluxes. Photo- 
metric redshifts were derived using EAZYv [34], including new custom templates 
designed to match the expected blue colors of early galaxies [35]. 


NOEMA Data and Analysis 


Observations with the JCMT and SCUBA2 bolometer array detected tentative 
850 um emission in the vicinity of CEERS-93316. [22] showed that such a 
detection is unlikely to come from a z > 10 object and explored different 
scenarios to explain the colors of this high-redshift galaxy candidate, finding 
that a z ~ 4.8 galaxy with some dust attenuation and strong emission lines 
could also explain the measured photometry. 

We observed dust continuum with band 3 of the Institute Radio Astronomie 
Millimétrique (IRAM) NOrthern Extended Millimeter Array (NOEMA) as a 
Director's Discretionary Time (DDT) program #D22AC001 (PI: S. Fujimoto). 
The observations were carried out between 2022 Nov 26 and 2022 Dec 08 in 
several visits with the C configuration using 10-11 antennas. The data were 


‘https: //ceers.github.io 
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processed in the standard manner with the pipeline using the latest version of 
the GILDAS software. We used CASA version 6.5 for the imaging [36]. 

The upper and lower side bands of the band 3 receiver were tuned to 
256.6 GHz and 271.1 GHz, resulting in the central wavelength of 1.134 mm. 
The observations of 3c345, 3c273, J1439--499, and 3c84 served as band-pass 
phase calibrators. Additional targets 1418--546 and J1439--499 were used for 
the phase and amplitude calibrations. We calibrated the absolute flux scale 
against MWC349, LkHal101, and 2010+723 whose flux densities are regularly 
monitored at NOEMA. The total integration time on-source was 6.5 hours. 

To maximize sensitivity, we used natural weighting for the imaging. The 
resulting 1 mm continuum map has synthesized beam FWHM of 0."91 x 0."67 
with lo sensitivities of 28 Jy beam! for the continuum. We do not identify 
any line features at the target position and thus produce the 1 mm continuum 
map from all channels. 

In Figure 5, we show the NOEMA 1.1-mm intensity contours overlaid on 
the NIRCam color image around CEERS-93316. We do not detect emission at 
the position of CEERS-93316, with a 3e upper limit of 0.084 mJy beam '!. 
This provides an upper limit of its dust-obscured star formation rate of ~20- 
200 Ms yr-! with a modified black body with dust temperature Ty = 40- 
70 K and a typical spectral index of 6g = 1.8, being consistent with the SFR 
estimate of 20+ 10 M /yr reported by [22]. This rules out the possibility that 
CEERS-93316 is a dusty starburst galaxy with SFR >> 200 Mo yr^!, while 
a moderately star-forming (SFR < 200 Mo yr~') dusty star-forming galaxy 
scenario remains. 

On the other hand, the NOEMA observations revealed relatively bright 
continuum emission from a nearby (~ 1."5) galaxy at 96 significance. We 
estimate 0.56+0.07 mJy with an optimized aperture. This is in good agreement 
with the previous report of SCUBA2 850 um emission of 0.64 + 0.26 mJy 
[22], implying that this neighboring galaxy is the main contributor to the 
aforementioned detection. 


Near-Infrared spectroscopy and data reduction 


JWST/NIRSpec multi-object spectroscopic observations (DDT program 
#2750; PI: P. Arrabal Haro) were executed on UTC 2023 March 24-25 using 
the PRISM/CLEAR combination of disperser and filter. The Micro Shutter 
Array [MSA; 37] was configured to observe 147 faint-object targets. The obser- 
vations consisted of 3 sequences of 3 nodded integrations, each with 9 to 10 
readout groups in NRSIRS2 readout mode, for a total combined exposure time 
of 18387 s. 

Three-shutter slitlets were employed, enabling a three-point nodding pat- 
tern to facilitate background subtraction. The prism disperser covers the 
wavelength range 0.6 to 5.3 um with varying spectral resolution R = A/AA = 
30 at A = 1.2,m to > 300 at A > 5m. The low resolution of the prism 
at bluer wavelengths aids the detection of faint UV continuum and the Lya 
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break, while the higher resolution at red wavelengths facilitates the detection 
of redshifted optical rest-frame emission lines. 

The MSA configurations included a total of 6 NIRCam-selected z = 8.5 
galaxy candidates as well as 141 other targets. We report results on 9 high- 
redshift candidates in Table 2 including the primary targets described above. 

The NIRSpec data processing followed the methodology employed for other 
CEERS NIRSpec prism observations [25]. We make use of the STScI Cali- 
bration Pipeline? version 1.8.5 and the Calibration Reference Data System 
(CRDS) mapping 1069. We use the calwebb_detector1 pipeline module to 
subtract the bias and the dark current, correct the 1/f noise, and generate 
count-rate maps (CRMs) from the uncalibrated images. The parameters of 
the jump step are modified for an improved correction of “snowball” events? 
associated with high-energy cosmic rays. 

The resulting CRMs are then processed with the calwebb_spec2 pipeline 
module which creates two-dimensional (2D) cutouts of the slitlets, performs 
the background subtraction making use of the 3-nod pattern, corrects the flat- 
fields, implements the wavelength and photometric calibrations and resamples 
the 2D spectra to correct the distortion of the spectral trace. The pathloss 
step accounts for the slit loss correction at this stage of the reduction process. 
Instead, when required for the analysis, we introduce slit loss corrections based 
on the morphology of the sources in the NIRCam bands and the location of 
the slitlet hosting them. 

The images of the three nods are combined at the calwebb_spec3 pipeline 
stage, making use of customized apertures for the extraction of the one- 
dimensional (1D) spectrum. The custom extraction apertures are visually 
defined for targets presenting high S/N at the continuum or emission lines 
when they are easily recognizable in the 2D spectra. For those cases where 
the targets are too faint for a robust visual identification, a 4-pixel extrac- 
tion aperture is defined around a spatial location estimated from the relative 
position of the target within its shutter, derived from the MSA configuration. 

The 2D and 1D spectra are simultaneously inspected with the Mosviz visu- 
alization tool* [38] to mask possible remaining hot pixels and other artifacts 
in the images, as well as the detector gap (when present). Data from the three 
consecutive exposure sequences are combined after masking image artifacts to 
produce the final 2D and 1D spectral products. 

The JWST pipeline uses an instrumental noise model to calculate flux 
errors for the extracted spectra. We test and rescale these flux errors for 
the effect of data correlation introduced by the pipeline following procedures 
described elsewhere [39]. 


?https://jwst- pipeline.readthedocs.io/en/latest /index.html 
3https://jwst-docs.stsci.edu/data-artifacts-and-features/snowballs-and-shower-artifacts 
^https://jdaviz.readthedocs.io/en/latest /mosviz/index.html 
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Redshift estimation 


In this paper, we focus on “high-priority” galaxies observed, which were 
selected based on either highly-likely high-redshift solutions, with one addi- 
tional source which showed a possible high-redshift solution but was more 
likely to be at lower redshift (MSA ID 2, discussed in [22]). For objects with 
multiple clear emission lines, we were easily able to identify the [O 11] and/or 
Ha emission lines. Spectroscopic redshifts are measured by comparing the 
observed emission lines with their theoretical vacuum transitions. A Gaussian 
profile is fitted for every spectral line to establish the transition wavelength. 
The Gaussian fluxes are used to weigh the mean redshifts in Table 2. In those 
observations, where the [O 111] 4960,5008 A doublet lines are merged, the pro- 
file fitting assumes the same gas kinematics and their amplitudes fixed by the 
theoretical emissivity ratio (2.98 see [40]) to isolate them. 

We also estimated the redshifts using the detected continuum with two 
methods. The first one consists in fitting the wavelength of the Lya break when 
it was observed. We do this following the method detailed in [39]. In sum- 
mary, we create a five-parameter model to describe an object/s spectrum. We 
first create the intrinsic spectrum based on the redshift, UV absolute magni- 
tude, and UV spectral slope. We then add IGM absorption by zeroing out the 
spectrum below 1215.67 À rest-frame wavelength, and add Lyo damping wing 
absorption by adding two additional parameters: the neutral hydrogen fraction 
(zur) and an ionized bubble radius (Rpubbie), following [41] and [18]. We derive 
posterior constraints on these five parameters using an IDL implementation of 
the EMCEE Python code (see [42] for details), restricting the spectrum to wave- 
lengths below 2500 À rest-frame for a given redshift. We performed a similar 
analysis using the MCMC sampler package from [43], finding broadly similar 
results, with most redshifts consistent to within Az ~ 0.1. The second method 
is full spectral fitting to the NIRSpec data, which detects the continuum for all 
our galaxies, using the EAZY [34] and a variety of spectral templates, including 
UV-bright and dusty galaxies. We remark that for Maisie's Galaxy, the spec- 
trum covers both the Lyman and Balmer breaks. All methods provide results 
consistent within Az ~ 0.1, consistent within our estimated uncertainties. 


Redshift results 


We derive emission-line redshifts for galaxies with MSA IDs 0 (CEERS-93316, 
[8]; via [O 11], [O 1m1], Ha), 2 ([22]; via [O rri], Ha), 28 (via [O 11], [O m]) and 355 
(via [O 111]), and Lya break redshifts for MSA IDs 1 (Maisie’s Galaxy; [12]), 28, 
64 and 355. Source ID 69 shows a weak break at ~1.3um, and this methodology 
does find a maximum likelihood of z — 9.744. However, the result has a large 
uncertainty, with a significant probability of the redshift being lower. We thus 
do not consider ID=69 confirmed. We present our redshift measurements in 
Table 2. 

For Maisie's Galaxy, we considered the alternative interpretation of a 
Balmer break at z — 3.114. Assuming that redshift and integrating the 
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observed spectrum in f, units through wavelength intervals 3145-3563 A and 
3751-4198 A [18] we measure signal-to-noise ratios of 0.0 and 9.9 below and 
above the break, and a 3c limit on the flux ratio > 2.6, which is considerably 
larger than would be predicted from stellar population synthesis models [18]. 
This supports the higher redshift interpretation of the break resulting from 
Lyo opacity of the IGM. 

We note that the derived spectroscopic redshifts are in broad agreement 
with the photometric redshift. For Maisie’s Galaxy, in particular, the pho- 
tometric redshift has been modified since the galaxy was discovered as data 
processing improved. The galaxy was originally thought to be at z ~ 14, but 
the photometric redshift was revised to z = 11.8*02 after improved astrometric 
registration recovered weak flux in F150W [7]. Improved NIRCam zero-point 
flux calibration revised this further to z = (WERE [12], while the latest photo- 
metric redshift from the photometry described in this paper (which is updated 
over that described in [12] to better match HST and JWST photometry) sug- 
gests z = 1108! oes consistent within the errors with our break-based redshift 
übs-—1L44* oe 

CEERS-93316 and the millimeter-detected galaxy CEERS-DSFG-1 (MSA 
ID 2) both have redshifts z = 4.91. We identify two more galaxies (MSA IDs 
2763 and 3149) at nearly the same redshift in the NIRSpec observations (Fig. 7, 
Table 2), supporting the hypothesis of a physical overdensity [20]. Moreover, 
the NOEMA-detected galaxy 1.5 arcsec from CEERS-93316 (Fig. 5) has a very 
similar photometric redshift and may well also be associated, but still requires 
spectroscopy. 


Emission lines 


The combined spectrum of “Maisie’s Galaxy” exhibits a particularly strong 
emission feature at ~ 4.6 um resembling a bright emission line. According 
to the redshift estimated from the Lya break, the only possible explanation 
for this line to be real is [O 11]. However, detailed analysis of individual nods 
revealed the presence of image defects that could be affecting this wavelength 
region. Hence, we do not consider this to be a real emission feature of the 
galaxy. Moreover, such a strong emission line being physical would necessarily 
imply four clear negative patterns above and below it in the 2D spectrum, a 
product of the nodded background subtraction and posterior combination (see 
the 2D spectra of MSA ID 28 and 355 in Figure 6). This is not the case for this 
hypothetical line, which suggests that it is indeed an artifact. All calculations 
performed with this spectrum were done after masking the spurious emission 
to avoid possible biases in the results it may cause. 

We also note that the observed strength of this putative [O 11] emission is 
not expected. This galaxy’s blue spectral slope implies a young, fairly low- 
metallicity stellar population (consistent with the low SED-fitting-derived dust 
attenuation). Such a stellar population should exhibit a fairly weak [O n1] line, 
with a high [Neimr]/[O 11] ratio [e.g., 26]. While the [Ne111|3869 line would 
be observed redder than [Or] where the noise is somewhat higher, should 
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the initially observed [O r1] flux be real, we would expect a significant detec- 
tion. However, the maximum signal-to-noise (SNR) within a few pixels of the 
expected position is ~1 after continuum subtraction. 

To explore what other constraints could be placed on emission lines in 
Maisie's Galaxy, we took the spectrum and subtracted the continuum emission 
as a rolling 40-pixel average. We then explored whether any pixels had a SNR 
of greater than two, finding no such pixels. We do find pixels with SNR just 
below 2 coincident with the expected observed wavelengths for O 1111660 and 
C 111]1909 at z = 11.383, which would be the spectroscopic redshift if the [O 1] 
emission were real. Given the low significance of these emission features, we do 
not consider these lines to be detected. However, if future deeper spectroscopy 
finds these to be robust, it would be consistent with our derived break redshift 
within the uncertainties. 

'The emission lines in CEERS-93316 were fitted manually from the spec- 
trum by first subtracting a local continuum and fitting Gaussian profiles to 
the visible emission lines. In the case of the blended but partially-resolved 
[O rui] 4960,5007 lines, a double Gaussian function was fitted. The Balmer lines 
(Ha and Hf) were corrected for underlying stellar absorption using the best- 
fitting stellar continuum from the Bagpipes fits (a « 296 correction in both 
cases; see below for a description of the Bagpipes fits). As described in the 
main text, these lines have large rest-frame equivalent widths (790+ 140À and 
1020 + 180A for [Om] and Ha, respectively), and thus significantly affect the 
long-wavelength NIRCam fluxes, mimicking the signature of a Lya break. 


Spectral energy distribution fitting 


We use three spectral energy distribution (SED) fitting codes to derive physical 
properties for our galaxies, Bagpipes [44], Cigale [45], and Synthesizer [46]. We 
apply several different methods to more thoroughly explore plausible values for 
the physical properties of these galaxies, as results obtained from SED fitting 
are known to be highly sensitive to the priors applied [47, 48], with different 
codes often returning a range of different values [49]. The photometric fluxes 
we fit are taken from [12], and we incorporate redshift constraints based on our 
new spectroscopic data. We describe in detail the approaches taken by each of 
these three codes separately below. 


Bagpipes 


Our Bagpipes fits make use of the 2016 updated version of the [50] stellar 
population models. Nebular emission lines and continuum are included in the 
model using results from the Cloudy photoionization code [51]. Dust atten- 
uation is included using a flexible-slope model, which also includes a Drude 
profile to model the 2175A dust bump and an additional factor multiplying the 
attenuation applied to star-forming regions (which we take to include all stel- 
lar emission arising from stars younger than 10 Myr and all resulting nebular 
emission) [52]. The details of our Bagpipes model assumptions are identical to 
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those used in [8, 53]. In particular, we use a constant star-formation history 
model, shown to produce stellar ages and stellar masses at the lower end of 
the plausible range for objects of interest. We apply a Gaussian prior to the 
redshift of each object based on our new spectroscopic data, with a standard 
deviation of 0.15 in each case. The median value for our Gaussian prior is set 
to z = 11.4 for Maisie's Galaxy and z = 4.9 for CEERS-93316. 


Cigale 


For this work, we use the last version of Cigale [54, 55], with updated emission 
line models computed from Cloudy [51]. For the dust, a Calzetti [56] extinction 
law is adopted for the dust attenuation of the stellar continuum. However, 
the nebular emission (continuum and lines) is attenuated with a screen model 
and a Small Magellanic Cloud (SMC) extinction curve [57]. No 2175À dust 
bump is added. We select a so-called “delayed + burst" SFH with SFR « 
t x exp(—t/r). Both the stellar and the gas metallicities are allowed to vary 
from low values Z = 0.0004 to Zo in the fit. As for Bagpipes, the mean value 
for our Gaussian prior is set to z — 11.4 for Maisie's Galaxy, and z — 4.9 
for CEERS-93316. For CEERS-93316, we use constraints from the equivalent 
widths measured on the NIRSpec spectrum of Ha, HG, and the two [O 11] 
lines. T'he derived equivalent widths estimated are in excellent agreement with 
0.7 < EW measurea/ EW CIGALE < 13. 


Synthesizer 


For the Synthesizer run, we explored star formation histories described by a 
delayed exponential with timescales ranging from 1 Myr to 100 Myr (i.e., from 
instantaneous to nearly constant star formation histories). We considered stel- 
lar populations models from [58], a [59] initial mass function, and metallicities 
smaller than the solar value. The nebular continuum and lines are included 
in the modeling as described in [46], with an attenuation law and differen- 
tial nebular-stellar dust effects as described in [60]. Redshifts are fixed to the 
values given in Table 2. For CEERS-93316, Ho and Hf equivalent widths 
were included in the fits, but all models provided lower values than those 
measured by a factor of ~ 4. For Maisie's Galaxy, a relatively strong degener- 
acy in the age was observed depending on the strength of the Balmer break, 
which is constrained by the spectroscopic data with limited SNR. Fitting data 
points bluewards of the Balmer break implied younger mass-weighted ages and 
smaller stellar masses (30410 Myr, log;; M./Mo = 8.3 X 0.1). Including NIR- 
Spec spectroscopic data points redwards of the Balmer break increased ages 
and masses (1202-30 Myr, logio M./Mo = 8.8 + 0.1). 


Data Availability. The JWST NIRSpec data are available from the Mikul- 
ski Archive for Space Telescopes (MAST; http://archive.stsci.edu), under 
program ID 2750. The CEERS JWST imaging data are available from MAST 
under program ID 1345. Reduced NIRCam data products from the CEERS 
team are available at https:/ /ceers.github.io. 
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Code Availability. JWST NIRSpec data were reduced using the 
JWST Pipeline (Version 1.8.5, reference mapping 1069; https:/ /github.com/ 
spacetelescope/jwst NIRSpec data inspection used the Mosviz visualization 
tool (https://jdaviz.readthedocs.io/en/latest /mosviz/index.html [38]. Photo- 
metric redshifts were measured using EAZY [34]. The Synthesizer code to 
perform stellar population synthesis modeling of galaxy spectro-photometric 
data is described in [46, 61]. The emission line profiles were fitted and their 
fluxes measured using LIME [62] (https://lime-stable.readthedocs.io). The 
Bagpipes code used for stellar population synthesis modeling is described in 
[44] and is available at (https://bagpipes.readthedocs.io). A detailed descrip- 
tion of the Cigale code can be found in [54, 55]. Cigale allows to fit data from 
the X-ray to the radio range and offers the possibility to add a contribution 
from an active galactic nucleus (AGN) and also to have constraints from non- 
photometric data (line fluxes, equivalent width, luminosities, or any physical 
parameter given in the input table). For this paper, we use a non-public ver- 
sion that fits photometric and spectroscopic data together. Cigale is available 
from (https://cigale.lam.fr/. 
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Extended Data 


r 


CEERS-93316 


Fig. 5 NOEMA 1.l-mm contours overlaid on the NIRCam 6” x 6” RGB color image 
(R: F444W, G: F356W, B: F200W) around CEERS-93316 marked with the white arrow. 
North is up and East is left. The magenta contours are drawn at lo interval from +20 to 
+90. We obtain the 3e upper limit of 0.084 mJy from CEERS-93316, while the continuum 
is robustly detected from a nearby source, which attributes to the tentative SCUBA2 flux 
(0.64 + 0.26 mJy [22]). 
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Fig. 6 2D and 1D NIRSpec spectra of five galaxies selected to have photometric redshifts 
z » 8, with spectral features and spectroscopic redshifts indicated, except for CEERS 13452 
where we are unable to determine a redshift with confidence. 
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Fig. 7 2D and 1D NIRSpec spectra of CEERS-93316 and three other galaxies at similar 
redshifts z ~ 4.9, highlighting a potential overdensity in the field at this redshift. 
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'Table 2: Summary of Spectroscopically Confirmed Galaxies 


MSA R.A. Dec F277W phot ens gpreak 
ID (J2000) (J2000) (AB) 68% C.L. 


0 214.914550 52.943023 26.5 16.0-16.6 4.912 + 0.001 — 4.89* 
1 214.943152 52.942442 27.9  10.7-11.5 — 11.44*009 11.38* 0-15 
2 214.909116 52.937205 26.2 3.4-13.6 4.910 + 0.003 — 4.88*0 
28 214.938642 52.011749 26.9 8.90-9.00 8.763 + 0.001 8.90'004  8.72*001 
64 214.922787 52.911529 28.3  10.4-11.6 = 10.1073 58 9.7310 


69 214.861602 52.904604 28.2 9.5-10.3 — — 
355 214.944766 52.931450 28.7 8.0-8.5 7.925 + 0.001 7.7143 43 7.887 


.01 
2763 214.927789 52.935859 25.1  4.87—5.44 4.902 + 0.001 = 4.87* 0.01 
3149 214.914917 52.943622 24.7 0.88-4.81 4.901 + 0.001 = 4.877 9 


Table 3: Summary of Confirmed Galaxy Physical Properties 


Software MSA ID Stellar Mass SFR Stellar Ay 
logio(M,/Mo) (Mo yr!) (mag) 
Bagpipes [44] 0 8.50.1 41 1.4 + 0.4 
1 8401 322 0.07* 0.09 
Cigale [54, 55] 0 9.5 + 0.2 8+ 2 1.8 + 0.3 
1 8.3 + 0.5 8+6 0.2 + 0.1 
Synthesizer [46] 0 8.8 + 0.1 60 + 20 2.3 + 0.2 
1 8.6 + 0.3 2-1 0.1 + 0.1 
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MSA ID 0 MSA ID 69° 


MSA ID 1 MSA ID 355 


MSA ID 64 MSA ID 28 


LE a 


Fig. 8 Montage of images of the galaxies with spectroscopic confirmation presented here. 
Each image shows a 3" x 3" image centered on the location of each galaxy, with ID labeled. 
The red-green-blue image for each galaxy corresponds to the JWST/NIRCam F444W, 
F356W, and F277W, respectively, rotated with North up and East to the left. The rectangles 
show the approximate location of the JWST/NIRSpec MSA slit positions. 
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